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Abstract 

We discuss the ways of extracting a low energy scale of an underlying theory using 
high energy scattering data. Within an exactly solvable model of quantum mechanics 
we analyze a technique based on introduction of nonperturbative power corrections 
accounting for asymptotically small terms and an alternative approach exploiting a 
modified running coupling constant of the model and nonperturbative continuation of 
evolution equations into an infrared region. Numerical estimates show that the latter 
is more efficient in approximating low-energy data of the model. 



PACS number(s): 11.15.Bt, 12.38. Cy 



1 Introduction. 



Considerable and steady improvement of experimental data has recently caused a renewal of 
interest in resummation of perturbation theory (PT) series. In the growing number of cases 
finite order theoretical predictions within PT have uncertainties comparable with experi- 
mental errors that urgently requires more accurate theoretical estimates. Nonperturbative 
(power-like) corrections to different processes within Wilson's OPE are widely used to match 
with improving of experimental data. For cases that have no simple formulation in terms 
of OPE the attempts to go beyond PT are now mostly based on using renormalons (for a 
concise up-to-date review see Because the expansion parameter - a running coupling 
constant as - is sufficiently large for moderate energies predictions differ strongly depending 
on a way one chooses to handle a strong coupling constant in the infrared region. 

The technique we further refer to as a standard one presently consists in resumming 
bubble chaines with principal value prescription for singularities in the Borel plane, e.g. 
0. Rich phenomenology can be developed on such a base though the real sensitivity 
of the approach to the infrared physics is unclear as well as an unambiguous disentangle 
of perturbative and nonperturbative (condensate) contributions 0]. Some other approaches 
use mostly the modified running of the coupling constant |[, specific recipes for scale 
setting 0, and optimization of perturbation theory |l^, or some modification of /5 
function to produce a smooth evolution at small momenta Initially there is no preference 



between these techniques because no exact results on the behaviour of PT in large orders 
or in IR domain are known. Some general properties of the quantum field theory to be 
respected (like analyticity) give no much help to distinguish between possibilities. Yet in 
phenomenological applications, the existing methods give different numerical results lying 
on the edge of experimental errors. The selection of a working frame will eventually be based 
on how well the particular technique fits experimental data. Under these circumstances it 
seems to be instructive to study some simple models |1^ where different methods used in 
QCD could be quantitatively checked for their advantages at least for gaining some intuition 
useful (or misleading as well) in more complicated situations. 

In the present paper we investigate two different approaches within a quantum mechanical 
model that mimics some general features of renormalons. 

2 The model. 

We consider the problem of potential scattering with 

V{r) = Vo6{r-ro) (1) 
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and limit ourselves to s-wave amplitudes. Such a potential can be considered as a kind of 
confining (not completely) one. We study a value of wave function at the origin (a free 
wave function is normalized to 1). The exact solution for scattering of the plain wave with 
momentum k reads 

^(fc) = ^{k, r = 0) = (^1 + y e^'^"" sin(A;ro)) . (2) 



where the mass of the particle in the hamiltonian is set to 1. To study scattering of wave 
packages with distributed momentum we consider an integral of the form 

PCX) 

$(A) = / ij{k)Wik,\)dk (3) 





where W{k,X) is a normalized weight function of a package depending on a set of parameters 
A, /o°^ W{k, X)dk = 1. It is more convenient to deal with a function F{X) 

^(A) = 1 + F(A) 

so that F vanishes if the scattering potential is switched off. Because of oscillating factors 
in eq. (0) integrals are not well suitable for the PT analysis (they are "Minkowskian" 
quantities). For m = r^^ > |Vo| there are no bound states in the potential (||) and we 
can carry out the Wick rotation because ipi^) is analytic in upper semi-plane (Im k > 0) 
that corresponds to the physical sheet in energy E k"^. In "Euclidean" region the exact 
expression for %lj{k) becomes 

^(g)= fl + ^(l-e-2«/-)) (4) 



where k = iq, q > 0. The last formula can be obtained by PT from a Born series for 
the standard Lippmann-Schwinger equation of potential scattering since we deal with a 
finite range potential. The Born series is ordered in Vq with general form of a nth term 
Vf^^n{q, ITT')- Bearing in mind the high energy analysis one can classify the behavior of Born 
series according the rate with which its terms vanish at large momenta. Then each term 
of the Born series for ip{q) ($n(g, m)) contains contributions vanishing as a power and as 
an exponent that is similar to the situation in QCD where PT contains two kinds of terms 
different with respect to their high energy behavior - logarithmic and power like. This 
different behavior allows one to easily separate terms accordingly. Such a structure of PT 
series is in fact a justification for us of using this model as a toy analog of QCD and one 
must always remember a toy character of the model. The parameter Vq then determines 
the scale at which the series of power vanishing terms becomes poorly convergent while the 
exponential terms have already died out completely. The parameter m determines the scale 
of exponentially supressed effects. In contrast to QCD where in the massless limit both 
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logatithmic and power suppressed terms are governed by a single dimensional parameter 
Aqcd; in our model one can change the relative weight of power and exponential terms 
varying parameters m and Vq independently or rather choosing the particular ratio Vo/m. 
Clearly at large values of this ratio exponetially small terms can hardly be detected at all 
and any modification of PT is going to be successful because corrections are really tiny. 
We however consider more interesting situation when this ratio is small enough and at 
some moderate energies both types of contributions are present with relatively essential 
weights. In such a situation one encounters a dilemma whether to keep them still separated 
therefore using the classification inspired by asymptotic considerations or to use more direct 
optimization technique though lifting the requirement of proper asymptotic behaviour but 
more efficient and precise at moderate energies. The numerical investigation of this problem 
is in fact the purpose of our paper. 

Note that we do not fix the sign of the parameter Vq and will study both attractive and 
repulsive interaction. 

At high energies {q ^ m) the expansion parameter is a running coupling constant a{q) = 
Vo/2g (trivial asymptotic freedom as in superrenormalizable theories) 

^lJ{q)=r'ia)+r''irn,Vo,q) 

where 

oo 

r'{a) = j:i-^r (5) 

n=0 

and ■ilj"'P{m, Vq, q) stands for exponentially suppressed "nonperturbative" terms 

rnm,\/o,g) = ^e-2- +••• (6) 

Within the present model we classify terms with respect to their behavior at large q: power 
like vanishing - PT, faster than any power - non-PT. With the standard renormalization 
group (RG) terminology we have 

P(«) = Q-g^ = (7) 

Resummation in eq. (in the spirit of RG) results in definition of a new (renormalization 
group improved) "running" coupling constant 

aas{q) 



2q + Vo 
with a (3 function 



naas) = -aas{l-a,s). (9) 
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This (3 function has an infrared fixed point that makes the running couphng constant aasiq) 
finite at small q in accordance with the explicit expression (|^). The use of this expansion 
parameter allows us to improve the perturbation theory and to sum up all "perturbative" 
power terms of the series 

r'{q) = l + aUq). (10) 

Now we turn to consideration of the wave package of a specific form given by the following 
weight function 

This weight function has a bump of the width v^<5 located at g ~ Q/2 so the above wave 
package can be considered as a "probe" of the scattering potential at the scale ~ 2/Q. 



3 Borel summation and renormalons. 

It is easy to see that our observable @ suffers from the renormalon. Substituting il){q) in 
eq. (H) by its asymptotic expansion (|ToD we obtain 



/■oo 

/ aUq)W{q,Q)dq. (11) 
Jo 



The quantity F°*(a) has a typical structure of QCD observable containing renormalon, i.e. 
it is an integral of some weight function multiplied by a running coupling constant over the 
interval that includes strong coupling domain. Note that in our model the use of the running 
coupling constant (H) in the integrand accounts for all perturbative corrections. Situation 
in QCD is much more comphcated and the representation of such a type is justified by the 



assumption of "naive nonabelianization" [13 



After integrating term by term in eq. ( [TTl) we get the series with factorially growing 
coefficients 

oo 

^"(«) = E(-l)"^!«"' a = K)/2<3. (12) 

n=l 

Properties of the series (0) depend crucially on the sign of a or Vq. Let us consider first 
repulsive potential Vq > 0. Then the alternating series (^2|) is Borel summable (in QCD it 
might correspond to an ultraviolet renormalon). The Borel image 

~ II 

F^\u) = --— (13) 

1 + M 

has a pole at u = — 1 and is a regular function on the positive semiaxis. So the Borel 
summation leads to the result (a = Vq/2Q) 

-l/a /^X 

Fb(«) = El - -1 (14) 

a \aj 
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where Ei(x) is the integral exponent [n 



Ei(x) = / 

J X 



OO 



e 
T 



Clearly the Borel resummation procedure in this case gives an unambiguous meaning to the 
series (^2]) and essensially improves on convergence of partial sums of the series (|12D . However 
numerical analysis shows that Fb{Q) does not approximate well the exact function F{Q) for 
intermediate Q m (Fig. 1). Indeed, aas{q) is the best expansion parameter (exponentially 
accurate) at large momenta but ip'^^{q) does not approximate well the function ipi^) at small q 
and eq. (0) tells us nothing about the parameter m that measures exponentially suppressed 
"nonperturbative" contributions. Is there more efficient way of extracting information on the 
parameter m? Alternative approach is to compute the function F{Q) within the modified 
perturbation theory for ip"'^{q) that can provide a sufficient accuracy even at very small q. 
For this purpose one has to choose a relevant expansion parameter. We write 

c^asiq) = a^\q) (l - ^««(g))" = ««(^) E i^X<^f{<lT. 



2q + fi 

where /x is a parameter that reflect some freedom in the choice of a scheme. Now we limit 
ourselves to only two terms of this expansion that is reasonable in PT region (though the 
precision may be improved on using more terms) and deflne new expansion parameter 



with a P function 



that is a PT transformation. The perturbation theory series for ^'{q) in a^^^ reads 

^-(g) = l + af (g)+0(a2(g)). (16) 

Taking the flrst order term in eq. (^) and flxing the parameter fi at an optimal value we 
flnd the function F^{Q) = C(^^\q)W{q, Q)dq to be very close to the exact function F(Q) 
up to very small Q (Fig. 1). Note that at very large Q the function Fb{Q) becomes closer 
to the exact result F{Q) than F^{Q) because Fb{Q) and F{Q) have the same asymptotic 
expansion by construction. On the other hand lifting this too strong condition of the same 
asymptotic behavior at inflnitely large q we flnd a function F^{Q) that approximates the 
exact function F{Q) uniformly for moderate Q. Because at the optimal /i the function ip^'^\q) 
is close to ip{q) for flnite q this approximation is universal in a sense that it works well for 
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various forms of scattering package. The optimal value /x°^*(m,Vo) can be extracted from 
experiment (in our model the exact solution plays the role of experimental data). It turns 
out to be very sensitive to the variation of m while the Borel resummed series is universal 
and gives the same prediction for theories with different m. So, in this case the modification 
of the running of the coupling constant is more efficient and flexible in determining the low 
energy structure of the model. Stress again that Borel resummed result in this case has no 
build-in mechanism to modify predictions in dependence on the parameter m while the value 
of optimal parameter /i directly and rather sensitively reflects the change of m. 

The case of attractive potential V{r) = —Vo6{r — ro), Vq > at flrst sight seems to be 
completely different. The running coupling constant 

has a singularity at g = Vo/2 > 0. The series ( [T^ ) now becomes 

oo 

F^'ia) = J2 ^'a"' " = '^o/2Q > (17) 

n=l 

and it is not Borel summable. This situation resembles the case of the infrared renormalon. 



The Borel image of the series ([T7|) has a pole on the positive semiaxis and the Borel procedure 
leads to an ill-deflned representation in case of attractive potential. Though this could be 
considered as a signal of the presence of nonperturbative contributions one should stress 
that the exact function (^) undergoes no qualitative change in the low energy domain (see 
also [|l^). Following the line of QCD renormalon technique we deflne the result of Borel 
summation by deforming the integration contour in the complex u plane. The result obtained 
in this way depends on the speciflc form of an integration contour while an appropriate 
nonperturbative part must cancel this dependence. In our model we use the principal value 
(PV) prescription to deflne the sum of the series (0) 

FM = PV r e-«/"^^ = ^Ei C-) - 1 (18) 
JO I — u a a \aj 



where Ei(x) is the integral exponent [|14 



/X pt 
-. 
-oo t 



a = Vq/2Q > 0. Within the renormalon technique one should search for the exact function 
F{Q) in the form 

F{Q) = Feia) + Ce"^/" + . . . (19) 

Here the flrst term has the same perturbative asymptotic expansion as the exact function 
F{Q), a constant C gives the leading exponentially suppressed correction and ellipsis stands 
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for "higher twist" contributions. The power of the exponent in ( ]T9[ ) is determined by the 
position of the pole of the Borel image. To find the value of C one has to use purely 
nonperturbative method or extract it from experiment. The form of first term, the value 
of the constant C, and high order corrections do depend on summation prescription while 
the whole sum should not by construction. So, in contrast to the previous case where the 
PT sum was uniquely determined by Borel prescription, in this case there is a parameter 
that can be adjusted to fit experimental data and therefore to measure the value of m. It 
happens that the parameter C(m, Vq) is quite sensitive to the variation of m so it can be 
considered "probe" of nonperturbative effects. 

The result of approximation for the exact function F{Q) with eq. (^) is given in Fig. 2. 

The approximation diverges strongly at small Q. The reason is the same as in the 
previous case of the Borel summable series. Namely, aas{q) is a poor expansion parameter 
at small momenta. Still at sufficiently large Q the accuracy is reasonable and adding of an 
adjusting term (constant C) improves on a precision of the resummed series at moderate 
energies. This is rather natural because an introduction of an additional free parameter (C 
in this case) always allows one to get better results. 

Though the running coupling aasiq) becomes singular at some point and can hardly be 
used for approximation of regular function ip{q) the introduction of a new term for Borel 
nonsummable series (extra degree of freedom) helps to detect the dependence of the constant 
C on m through experimental data. Note that this method gives no hint on how to construct 
"higher twist" contributions in this case because there is a single singularity of the Borel 
image on the positive semiaxes allowing only one additional degree of freedom. 

Stress once again however that even the singular function aas^q) still accumulates all 
PT terms exactly as in the previous case. Though making results numerically better, in our 
model the assumption that a singularity of the Borel image determines the index of vanishing 
of the leading nonperturbative term is incorrect. Indeed eq. (p!9[) fails to detect the leading 
non-PT asymptotics of the exact function F{Q) which reads 



while the Borel based guess is exp{—l/a), a = Vq/2Q. In addition, the parameterization (|T9| ) 
is not universal i.e. one gets essentially different values of C for different scattering packages 
and it is not clear how this variety of parameters should be used to estimate the infrared 
parameter m of the model. 

Turning to an alternative approach and introducing a running coupling constant of the 
form (^) with an appropriate fi we obtain a uniform approximation of the function F{Q) 
practically for all Q (Fig. 2). Because in our toy model experimental data (exact solu- 
tion) is available for all Q with arbitrary accuracy we do not limit the applicability of the 
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methods in question to the asymptotic regime only but try them at all possible Q. Thus, 
the approach based on introduction of a single parameter as in the previous case produces 
more efficient fitting for experimental data in broader range of energies. The accuracy of 
the approximation and the actual value of the parameter /i (and C as well) depend on the 
range where data have to be ffited. The renormalon approach cannot be used for small Q by 
construction of approximants while the modified running technique is applicable till rather 
small energies. The optimal value of the parameter /x depends on the scattering packages 
rather weakly because this technique essentially approximates the exact function point- 
wise. This dependence is also expected because it reflects the different choice of optimal 
scheme for different observables. 

In QCD this prescription would correspond to the use of (probably mass dependent) RG 
equation for strong coupling constant with infrared regular solution |TT], |T^. We should 
note however that in this way we determine only the scale of nonperturbative effects while 
the exact form of F{Q) can be found only via real nonperturbative calculations. 

4 Conclusion. 

The analysis of two possible ways of extracting information on low energy domain of a 
quantum mechanical model shows that approach based on optimization of the PT through 
introduction of a flexible expansion parameter using the freedom of choice of the scheme is 
more efficient for moderate energies than direct Borel resummation technique. In our model 
also the singularity of Borel image for attractive potential does not correspond to leading 
non-PT asymptotics of exact function that is one of main reasons for using the renormalon 
technique in phenomenological applications of QCD. Though obtained in a toy model, these 
observations may serve as a ground for using a modified running of the coupling constant of 
QCD in the infrared domain for phenomenological applications instead (or in addition to) 
the standard renormalon technique. 
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Figure Captions 

Fig. 1. Numerical results for repulsive potential Vq = 1, m = 3. 
Function F^^{Q) j F{Q)^ the optimal value /x = 5 (line a). 
Function Fb{Q)/F{Q) (line b). 

Function F^^\Q)/F{Q) where F^^\Q) is the sum of the first three terms of the asymptotic 
expansion ( |T2|) (line c). 

Fig. 2. Numerical results for atractive potential Vq = 1, m = 3. 
Function F^^{Q) / F{Q) for the optimal value = 2.7 (line a). 
Function {Fb{Q) + Ce~^/'') / F{Q) for an optimal value C = -0.06 (line b). 
Function Fb{Q)/F{Q) (line c). 



11 




0123456789 10 

g/2 

Fig. 1 



12 




Fig. 2 
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